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Refined astrometry measurements allow us to detect large-scale deviations from isotropy through 
real-time observations of changes in the angular separation between sources at cosmic distances. 
This "cosmic parallax" effect is a powerful consistency test of FRW metric and may set independent 
constraints on cosmic anisotropy. We apply this novel general test to LTB cosmologies with off- 
center observers and show that future satellite missions such as Gaia might achieve accuracies that 
would put limits on the off-center distance which are competitive with CMB dipole constraints. 



Introduction. The standard model of cosmology rests 
on two main assumptions: general relativity and a homo- 
geneous and isotropic metric, the Friedmann-Robertson- 
Walker metric (henceforth FRW). While general relativ- 
ity has been tested with great precision at least in lab- 
oratory and in the solar system, the issue of large-scale 
deviations from homogeneity and isotropy is much less 
settled. There is by now an abundant literature on tests 
of the FRW metric, and on alternative models invoked to 
explain the accelerated expansion by the effect of strong, 
large-scale deviations from homogeneity (see [T] for a re- 
view). Several such models adopt as an alternative the 
Lemaitre-Tolman-Bondi metric (henceforth LTB), that 
is a model with a spherically symmetric distribution of 

matter (see for instance [2]|3l[I]|5]|6j[7]|8]). ma ^ n m0 ~ 
tivation for this is the fact that the distance-dependent 
expansion rate can explain the supernovae la excess dim- 
ming without a dark energy field. 

LTB universes appear anisotropic to any observer ex- 
cept the central one. In every anisotropic expansion 
the angular separation between any two sources varies 
in time, thereby inducing a cosmic parallax (CP) ef- 
fect. This is totally analogous to the classical stellar 
parallax, except here the parallax is induced by a dif- 
ferential cosmic expansion rather than by the observer's 
own movement. Together with the Sandage effect of 
velocity shift z j9j [10], CP belongs to the new realm 
of real-time cosmology, a direct way of testing our uni- 
verse based on cosmological observations spaced by sev- 
eral years [TT1 IT2l IT3l PHI IT5] . 

One can expect on dimensional grounds that this dif- 
ferential cosmic expansion generates after an observation 
time lag of At = 10 years an effect of order HoAt = 
10~ 9 h; therefore two sources separated by 1 rad today 
will show a parallax of 10~ 9 h rad ps 200ft. /ias in ten 
years, which is well above the accuracy goal of 10 /xas 
set by Gaia (for V<15) [16] and other planned mis- 
sions like SIM H2, JASMINE [18] and VSOP-2 [19]. 
Adopting Gaia specifications, we show that 10 6 quasars 
might be enough to constrain the off-center distance ro 
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Figure 1: Overview, notation and conventions. For clarity 
purposes we assumed here that the points C, 0, 01, 61, ffl2, &2 
all lie on the same plane. Comoving coordinates r and ro 
correspond to physical coordinates X and Xq. 

to the 10 Megaparsec scale, similar to or better than 
any other proposed test of the Copernican principle but 
without the degeneracy with our peculiar velocity that 
afflicts the constraints from the cosmic microwave back- 
ground [5] 8, 20] . In this paper we first derive an estimate 
of CP by intuitive arguments based on a FRW descrip- 
tion and then test it with a full numerical integration of 
light-ray geodesies in LTB models. 

Estimating the cosmic parallax. Figure [l] depicts the 
overall scheme describing a possible time-variation of 
the angular position of a pair of sources that expand 
anisotropically with respect to the observer. We label 
the two sources a and b, and the two observation times 
1 and 2. In what follows, we will refer to (t, r, 9, tp) as 
the comoving coordinates with origin on the center of a 
spherically symmetric model. Peculiar velocities apart, 
the symmetry of such a model forces objects to expand 
radially outwards, keeping r, 8 and <f> constant. 

Let us assume now an expansion in a flat FRW space 
from a "center" C observed by an off-center observer O 
at a distance Xq from C. Since we are assuming FRW 
it is clear that any point in space could be considered a 
"center" of expansion: it is only when we will consider a 
LTB universe that the center acquires an absolute mean- 
ing. The relation between the observer line-of-sight angle 
£ and the coordinates of a source located at a radial dis- 
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tance X and angle in the C-frame is 

X cos - X a 



cos£ = 



(1) 



{X 2 + X 2 -2X Xcos0y/ 21 

where all angles are measured with respect to the CO axis 
and all distances in this section are to be understood as 
physical distances. 

We consider first two sources at location a% , b\ on the 
same plane that includes the CO axis with an angu- 
lar separation 71 as seen from O, both at distance X 
from C (throughout this letter we shall always assume 
for simplicity that both sources share the same <j> coor- 
dinate). After some time At, the sources move to po- 
sitions 02,62 and the distances X and Xq will have in- 
creased by A t X and A t Xo respectively. If for a moment 
we allow ourselves the liberty of assigning to the scale 
factor a(t) and the H function a spatial dependence, a 
time-variation of 7 is induced. The variation A t j is the 
anticipated cosmic parallax effect and can be easily es- 
timated if we suppose that the Hubble law is just gen- 
eralized to A t X = XH(t ,X)At = XH x At. General- 
izing to sources on different shells separated by a small 
AX = Xfy — X a (not to be mistaken with the time in- 
terval A t X) and AO = O^i — 6 a \, after straightforward 
geometry we arrive at 

A t7 = s At(H ohs - H x ) (cos 6 AO + sin AX/X) 
+ H x dH x /dX sin 6 AX + 0{s 2 ) . 

where if b s = H(to,ro), s = Xo/X <C 1 (at this or- 
der we can neglect the difference between the observed 
angle £ and 0). We can also convert the above inter- 
vals AX into the redshift interval Az by using the rela- 
tions r = J dz/H(z) and X = J a(t ,r)dr, which com- 
bine to AX = a(to, X)Az/H(z) ~ Az/H(z) (we impose 
a(t ,X ) = 1), where H(z) = H(t(z),X). 

In a FRW metric, H does not depend on r and the par- 
allax vanishes. On the other hand, any deviation from 
FRW entails such spatial dependence and the emergence 
of cosmic parallax, except possibly for special observers 
(such as the center of LTB). A constraint on A t 7 is there- 
fore a constraint on cosmic anisotropy. 

Rigorously, one actually needs to perform a full inte- 
gration of light-ray geodesies in the new metric. Never- 
theless, we shall assume for a moment that for an order 
of magnitude estimate we can simply replace H with its 
space-dependent counterpart given by LTB models. In 
order for a LTB cosmology to have any substantial effect 
(e.g., explaining the SNIa Hubble diagram) it is reason- 
able to assume a difference between the local H ^ s and 
the distant Hx of order 7? b s IS]. More precisely, putting 
H b s — Hx = H f, s Ah then one has that the overall effect 
A t 7 for on-shell sources is of order scos 0H o b s AhAtAO. 
That is, as anticipated, for At = 10 years we expect 
a parallax of 200s cos OAhAO //as for sources separated 
by AO. Similarly, for radial source pairs, neglecting 
the dHx/dX term (which is valid except close to the 



edge of the LTB void), one has ssm0H o b s AhAtAzjz « 
200s sin# AhAz/z /ia,s (assuming X ss zH(z)~ 1 ). 

Let us finally consider the main expected source of 
noise, the intrinsic peculiar velocities of the sources. The 
variation in angular separation for sources at angular di- 
ameter distance Da (measured by the observer) and pe- 
can be estimated as 



culiar velocity v pcc 



A(7pe 



500km/sy \lGpcy \ 10 years 



At 



lias. (3) 



This velocity field noise is therefore typically smaller than 
the experimental uncertainty (especially for large dis- 
tances) and again will be averaged out for many sources. 
Notice that the observer's own peculiar velocity produces 
a systematic offset sinusoidal signal A t 7 pcCj o of the same 
amplitude as A t 7 pec that has to be subtracted from the 
observations: we discuss this further below. 

Geodesic Equations. These very suggestive but simplis- 
tic calculations need confirmation from an exact treat- 
ment where the full relativistic propagation of light rays 
is taken into account. We will thus consider in what fol- 
lows two particular LTB models capable of fitting the 
observed SNIa Hubble diagram and the CMB first peak 
position and compatible with the COBE results of the 
CMB dipole anisotropy, as long as the observer is within 
around 15 Mpc from the center [5]. Moreover, both mod- 
els have void sizes which are small enough (z ~ 0.3) not 
to be ruled out due to distortions of the CMB blackbody 
radiation spectrum [20 . 

The LTB metric can be written as (primes and dots 
refer to partial space and time derivatives, respectively): 



ds 2 = -dt 2 



-dr 2 + R 2 (t,r)dn 2 



(4) 



[R'(t,r)Y 
l + /3(r) 

where (3(r) can be loosely thought as position dependent 
spatial curvature term. Two distinct Hubble parameters 
corresponding to the radial and perpendicular directions 
of expansion are defined as H\\ = R' /R' and H± = R/R 
(in a FRW metric R = ra(t) and H\\ = H±). This class of 
models exhibits implicit analytic solutions of the Einstein 
equations in the case of a matter-dominated universe, to 
wit (in terms of a parameter rf) 
ex 

R = (cosh 77 - 1)— + R h 



\ff3t = (sinh 77 



'2/? 
V) a/(2(3) - 



cosh rj + v D sinh r\ , (5) 



R\ ss sinh rj + \f~D (cosh 77 — 1) 



(6) 



where D — (a + (3R\ SS ) /((3Ri ss ), and a, (3 and i?i ss are 
all functions of r. In fact, R\ ss (r) stands for R(0,r) 
and we will choose t = to correspond to the time of 
last scattering, while a(r) is an arbitrary function and 
/3(r) is assumed to be positive. Due to the axial symme- 
try and the fact that photons follow a path which pre- 
serves the 4- velocity identity u a u a = 0, the four second- 
order geodesic equations for (t, r, 0, (j>) can be written as 
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five first-order ones. We will choose as variables the 
center-based coordinates t, r, 6, p = dr/dX and the red- 
shift z, where A is the affine parameter of the geodesies. 
We shall refer also to the conserved angular momentum 
J = R 2 d8/d\ = const = Jq. For a particular source, the 
angle £ is the coordinate equivalent to 9 for the observer, 
and in particular £o is the coordinate £ of a photon that 
arrives at the observer at the time of observation to- Ob- 
viously this coincides with the measured position in the 
sky of such a source at i . In terms of these variables, and 
defining A such that it(A) < 0, the autonomous system 
governing the geodesies is written as 



dt 
dX 

d z 
dX 

dp 
dX 




J 



(7) 



Following [5], the angle £ along a geodesic is given by 
cos£ = —R'(t,r)p/(uy/l + (3{r)), from which we ob- 
tain po = + P(r ) cos(£ ) / R'(t , r ) and J = J = 
R(to, To) sin(£o)- Therefore, our system is completely de- 
fined by the initial conditions t , r , 8 = 0, z = and 
£o- The first two define the instant of measurement and 
the offset between observer and center, while £o stands 
for the direction of incidence of the photons. By inte- 
grating the geodesic equations for two sources located at 
(z a i , Zbi , £<ji , £&i ) after a time interval At the CP will be 

A t 7 = 72 - 7l = (£a2 - 62) - (£al - &>l)- 

The models of Ref. [5] are characterized by a smooth 
transition between an inner void and an outer region with 
higher matter density and are described by the functions 
a(r) and (J(r) (Eqs. (28) and (29) in (5|), themselves 
carrying a total of 4 free parameters, one of which the 
value of the Hubble constant at the outer region, 

set at 51 km/ (s Mpc) . Following [5] we dub them Model 

I and II; the main difference between them is that Model 

II features a sharper transition from the void. How- 
ever transition width is not expected to be an impor- 
tant factor in CP since most quasars are outside the void 
and the most relevant quantity is the difference between 
the inner and outer values of H. In both cases we set 
the off-center (physical) distance to 15 Mpc, which is 
the upper limit allowed by CMB dipole distortions [5], 
and this corresponds to s ~ 4.4 1CP 3 for a source at 
2=1. It can be shown that in a step-like LTB void 

, A vo / X - 



H^(l- 



model, the Hx in Q is given by Hu Q 
X vo /X) a fig* whereas H ohs a H^ Q . 

In Figure [2] we plot A t j for three sources at z = 1, 
for models I and II as well as for the FRW-like esti- 




Figure 2: Atj for three sources at the same shell, at z = 1, for 
both models I (full lines) and II (dashed), and the FRW-like 
estimate (dotted). As expected, At 7 is linear in A9. 



mate. One can see that the results do not depend sen- 
sitively on the details of the shell transition and that 
in both cases the FRW-like estimate gives a reasonable 
idea of the true LTB behavior. We conclude that ([2j 
is a valid approximation. Numerically, we find that a 
convenient estimate of the parallax is given by A t 7 = 
100 sAh(At/W years) x (A6>, Az/z)^as. For the LTB 
models above Ah « 0.12 - 0.14. 

Integrating Q yields also another interesting observ- 
able, z = Atz/At, coupled to the CP: together they form 
a new set of real-time cosmic observables. Ref. [TJ 3 calcu- 
lated i for an observer at the center of a LTB model. In 
the limit of small s our numerical results reveal that z for 
off-center observers show very small angular dependence. 

Cosmic parallax with Gaia. A realistic possibility of 
observing the CP is offered by the forthcoming Gaia 
mission. Gaia will produce in five years a full-sky map 
of roughly 500, 000 quasars; by making of order 100 re- 
peated measurements over the five years mission, Gaia 
will hopefully achieve a positional error p between 10 to 
200 /ias (for quasars with magnitude V = 15 to 20). To 
compare our observations to Gaia we need to evaluate the 
average A t 7 with At = 5 years and N sources. The final 
Gaia error p is obtained by best-fitting 2N independent 
coordinates from N 2 /2 angular separation measures; the 
average positional error on the entire sky will scale there- 
fore as (2iV)" 1 / 2 . Over one hemisphere we can therefore 
estimate that the error scales as p/vN . Since the av- 
erage angular separation of random points on a sphere 
is 7r/2, the average of A t 7 can be estimated simply as 
Atj(0 = 7r/2). We find numerically A t 7(7r/2) = 10 s /zas, 
with little dependence on Az. Therefore Gaia can see the 
parallax if p/y/N < lOs/xas. For s = 4.4- 10~ 3 (i.e. the 
current CMB limit) and p = 30 ^as we need N > 450, 000 
sources: this shows that Gaia can constrain the cosmic 
anisotropy to CMB levels. An enhanced Gaia mission 
with At = 10 years (or two missions 5 years apart), 
p = 10 /ias and N — 10 6 would give s < 5 • 10~ 4 , i.e. 
ro < 1 Mpc if we assume the sources are at 2 Gpc. Two 
caveats of this preliminary estimate are however to be 
kept in mind. First, to map the angular change, Gaia's 
observing strategy should be redesigned in order to maxi- 
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mize the time interval between quasar observations. Sec- 
ond, being planned to monitor the local matter distribu- 
tion, Gaia's errors estimates are based on using a fraction 
of the quasars as a reference frame. Observing the quasar 
CP would require a different statistical analysis. Both ef- 
fects are likely to increase the expected final error. 

Two local effects induce spurious parallaxes: one (of 
the order of 0.1/ias yr~ v ) is induced by our own pecu- 
liar velocity and the other (of the order of 4/xasyr _1 [2T] 
by a changing aberration. Both produce a dipolar sig- 
nal, just like a LTB: however, the peculiar velocity par- 
allax decreases monotonically with the angular diameter 
distance, while the aberration change is independent of 
distance [21]. In contrast, the LTB signal has a character- 
istic non-trivial dependence on redshift: for the models 
investigated here it is vanishingly small inside the void, 
large near the edge, decreasing at large distances. It is 
therefore possible in principle to subtract the cosmic sig- 
nal from the local one, for instance estimating the local 
effects from sources inside the void, including Milky Way 
stars. A detailed calculation needs a careful simulation 
of experimental settings (including possibly effects like 
source photocenter jitter and relativistic light deflection 
by solar system bodies) which is outside the scope of this 
paper. Moreover, more general anisotropic models will 
not produce a simple dipole. 

Conclusions. Planned space-based astrometric mis- 
sions aim at accuracies of the order of few microarcsec- 
onds. In this paper we have shown that the cosmic paral- 
lax of distant sources in an LTB model might be observ- 
able employing the same missions. Similar considerations 
would apply to all other anisotropic cosmological mod- 
els as e.g. the Bianchi models. A positive detection of 
large-scale CP would disprove therefore one of the basic 
tenets of modern cosmology, isotropy. We have shown 
that for a typical LTB model designed to explain the 
supernovae la Hubble diagram a enhanced Gaia experi- 
ment could constrain the anisotropy parameter s to less 
than 10~ 3 , corresponding to the Megaparsec scale, much 
better than current CMB dipole limits [5]. Moreover, 
this test may probe a different range of scales depending 
on the quasar redshift distribution and, contrary to the 
CMB limits, the CP method cannot be completely under- 
mined by the observer's peculiar motion and is limited 
only by source statistics instead of by the cosmic vari- 
ance. We anticipate however that the major source of 
systematics would be the subtraction of the aberration 
change parallax. Real-time cosmology tests directly cos- 
mic kinematics by observing changes in source positions 
and velocities. We have shown that the cosmic parallax, 
along with the velocity shift effect i, can fully reconstruct 
the 3D cosmic flow of distant sources. 
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